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ABSTRACT 



Thermonuclear runaways driven by accretion into degenerate white dwarf 
cores are the source that power classical nova outbursts. In this paper, we iden- 
tify the dominant nuclear paths involved in the synthesis of intermediate-mass 
elements, from Si to Ca, during such violent events. New evolutionary sequences 
of 1.35 Mq ONe novae have been computed, using updated nuclear reaction 
rates. The main nuclear activity in this region is powered by the leakage from 
the NeNa-MgAl region, where the activity is confined during the early stages of 
the explosion. We discuss the critical role played by ^°P(p,7) in the synthesis of 
nuclear species beyond sulfur and point out the large uncertainty that affects its 
rate, which has dramatic consequences for studies of nova nucleosynthesis in the 
Si-Ca mass region. 

Subject headings: novae, cataclysmic variables — nuclear reactions, nucleosyn- 
thesis, abundances 
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1. Introduction 

Classical novae inject nuclear processed material into the interstellar medium, and there- 
fore may play a role in the Galactic nucleosynthesis. Isotopes likely attributed to nova 
outbursts include ^^C, ^^N and ^''O (see Jose & Hernanz 1998, and references therein). Cur- 
rent hydrodynamic models of classical novae show that explosions on carbon-oxygen white 
dwarfs (hereafter, CO novae) are dominated by the synthesis of CNO-group nuclei, largely 
overproduced in the ejecta with respect to solar abundances (Starrfield et al. 1997; Kovetz 
& Prialnik 1997; Jose & Hernanz 1998). On the other hand, explosions on oxygen- neon 
white dwarfs (hereafter, ONe novae), more massive than CO ones, provide large enrichments 
not only in CNO nuclei, but also in neon, sodium, aluminum, or other intermediate-mass 
elements (Starrfield et al. 1998; Jose & Hernanz 1998; Jose et al. 1999). Hence, the com- 
position of the ejecta is a key signature of the nature of the underlying white dwarf (CO 
or ONe) and provides also valuable information on the explosion (i.e., peak temperature, 
characteristic timescale...). In particular, the presence of isotopes of the silicon to calcium 
mass region (hereafter, Si-Ca), overproduced with respect to solar values, in the ejecta of 
a classical nova, may reveal the presence of an underlying massive white dwarf (the ones 
reaching the higher peak temperatures during the explosion) since proton captures on nuclei 
above silicon require temperatures in excess of 3 x 10^ K to overcome their large Coulomb 
barriers. Indeed, abundance determinations from some observed novae, such as V1370 Aql 
1982 (Snijders et al. 1987), QU Vul 1984 (Andrea et al. 1994), V2214 Oph 1988 (Andrea et 
al. 1994), or V838 Her 1991 (Vanlandingham et al. 1996, 1997), show the presence of Si-Ca 
nuclei in the ejecta. 

Moreover, the nova contribution to the Galactic content of other species such as ^Li 
(Starrfield et al. 1978, Hernanz et al. 1996) or ^^Al (Pohtano et al. 1995, Jose et al. 1997, 
Starrfield et al. 1998, Jose et al. 1999) may also be significant, although its real extent is 
still a matter of debate (see Romano et al. 1999, for details on the contribution of novae to 
match the ^Li content in realistic calculations of Galactic Chemical evolution). The synthesis 
of such medium- or long-lived species has, in turn, implications for the high-energy emission 
during nova outbursts: several gamma-ray signatures have been predicted for nova outbursts 
(Clayton & Hoyle 1974; Leising & Clayton 1987; Gomez-Gomar et al. 1998, Hernanz et al. 
1999), including a 478 keV line from ''Li (CO novae), a 1275 keV line from ^^Na (ONe novae), 
and an intense 511 keV line and the corresponding continuum at energies below (for both CO 
and ONe nova types). In the Si-Ca mass region, another gamma-ray signal, coming from 
the P'^-decay of the medium lived ^^"*C1 (isomeric state), has also been proposed for 511 
keV (Leising & Clayton 1987) and 2.128 MeV (Coc et al. 2000a) line emission. However, 
no nucleosynthesis calculation with realistic nuclear reaction rates has, up to now, been 
attempted to evaluate the chances for the potential detection of this particular gamma-ray 



-3- 



signal from novae. 

Furthermore, the recent discovery of a number of presolar grains hkely formed in the 
ejecta of classical novae (Amari et al., 2001) also points out the importance of a good 
determination of the abundance pattern, in particular in the Si-Ca mass region: together 
with a number of constraints from specific isotopic ratios (i.e., low -"^^C/^^C and ^^N/^^N 
values, high ^^Al/^^Al), the basic feature that points towards a nova origin relics on the large 
excesses in ^°Si/^^Si and close-to-solar ^^Si/^^Si measured in such grains. The determination 
of other isotopic ratios in this mass range, such as ^^S/^^S, may also help to identify future 
nova grain candidates. 

These aspects point out the interest of a more accurate determination of the relevant 
synthesis mechanisms in the Si-Ca mass region. Such an analysis has been only scarcely 
attempted (see recent parametrized nucleosynthesis calculations by Lefevbre et al. 1997, 
Wanajo et al. 1999, and Iliadis et al. 1999, as well as hydrodynamic calculations by Starrfield 
et al. 1998), partially because of the large number of isotopes of interest in this mass 
region, around 40 species, and the hundred nuclear reactions involved. Accordingly, it is 
the purpose of this paper to determine the dominant nuclear paths involved in the synthesis 
of intermediate-mass elements, from Si to Ca, during classical nova outbursts, identifying 
the key reactions and investigating their associated nuclear uncertainties, that may have a 
dramatic influence in the overall process. In Section 2, we outline the basic nuclear paths 
at different stages of the nova outburst for a particular model. Details on nuclear physics 
aspects are given in Section 3. Constraints on the production of Si-Ca group nuclei, together 
with the most relevant conclusions of this paper are summarized in Section 4. 



2. Main nuclear paths in the synthesis of Si-Ca nuclei in classical novae 

The high temperatures required to initiate a noticeable nuclear activity in the Si-Ca 
group nuclei restrict the domain of interest to explosions in very massive white dwarfs (i.e., 
close to the Chandrasekhar mass). In this section, we will focus on the main nuclear paths 
involved in the synthesis of several isotopes within this mass range, through a detailed 
analysis of a 1.35 ONe white dwarf, which accretes solar-hke matter at a rate of M = 
2 X lO"^'^ Mg.yr"^ , assuming a 50% degree of mixing with the outermost shells of the 
white dwarf core. Initial abundances of the Si-Ca group nuclei adopted in this Model are 
listed in Table 1 (abundances up to Si are given in Jose et al. 1999, Table 1). In fact, the 
specific prescription adopted for the initial abundance pattern is crucial for nucleosynthesis 
calculations. Whereas a deep analysis of this effect is out of the scope of this paper, and will 
be presented elsewhere, we just point out that it introduces some uncertainty in all current 
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nucleosynthesis calculations during nova outbursts. 

The main source for the production of the Si-Ca isotopes (see Fig. 1) comes from the 
processing of the initial Ne-Al isotopes from the white dwarf leaking out of the NeNa and 
MgAl mass regions, mainly through ^^Si. We do not discuss here the nuclear processing 
within these regions but refer the reader to a previous paper (Jose et al. 1999). Hence we 
start our discussion with the last reactions linking the MgAl to SiP mass regions. The new 
evolutionary sequences of nova outbursts have been calculated by means of the code SHIVA 
(see Jose & Hernanz 1998 for details), a one-dimcnsional implicit hydrodynamical code that 
follows the course of the explosion from the onset of accretion up to the ejection stages. 
Calculations have been performed with a nuclear reaction network updated for the region of 
interest (see section 3) . Snapshots of the evolution of the most relevant elements are shown 
in Figures 2 to 6. 

At the early stages of the explosion, nuclear activity above CNO cycle is mostly confined 
in the NeNa and MgAl region with only marginal activity in the Si-Ca region. Soon after 
the thermonuclear runaway takes place, when the temperature at the burning shell reaches 
Tf,s — Sx 10^ K (Figs. 2 to 6, first panel), the main nuclear path affecting Si to Ca evolution is 
dominated by a series of reactions, beginning with proton captures onto ^^^Al (ground state), 
2^5Al(p,7)27Si, followed by '^'^Si{(3+f'^Al{p,'yf^Si. As a result, ^^Si is slightly enhanced at the 
burning shell, since its /?+-dccay channel is slower than ^^^Al(p,7). Marginal nuclear activity 
at the burning shell involves other proton-capture reactions on ^"'^Al and ^Sm^j (isomeric 
state), as well as onto ^^Si, which account for a moderate increase in ^^Si and ^^P. No 
relevant nuclear activity beyond phosphorus is found. 

Evolution at T^s = 10^ K (Figs. 2 to 6, second panel) is still dominated by the chain 
263Al(p, 7)27Si(/?+)27Al(p,7)28Si, plus the contribution of ^^A\{p,-ff^Si{(3+f^'^A\{p,-ff'^Si, 
and ^'^Si(p,7)^^P(/9"'")^^Si. Some marginal activity involves also '^"'^P(p,7)'^^S, and ■^'^S(p,7)^^Cl. 
The most remarkable feature at this stage is an important increase in both ^^'^^Si. The 
outermost shells are dominated by some marginal decays (the lower temperatures prevent 
proton captures to proceed efficiently), such as '^^Si{P+f^A\, '^^Si{P+f^"'A\, and =^°P(/3+)=^°Si, 
followed by 29p(/3+)29Si, and 34"»Cl(/3+)=^^S. 

Important changes take place as soon as T^g reaches 2 x 10^ K (Figs. 2 to 6, third panel). 
The nuclear activity is now driven by ^^^Al(p,7)2^Si(p,7)2®P (faster than ^'^Si{(3+y'^ Al), 
and by ^'^ A\{p,^y^Si, followed by 25Al(p,7)26Si(/3+)26"^Al(p, ^f^SiiP+y^Al, 28p(/3+)28Si, 
26Si(p,7)27p(/3+)27si, 28p(p^^)29g^ 28,29gj(p^^)29,30p_ rpj^g outermost shclls are dominated 
by 27Si(/5+)27Al, 26Si(,5+)26-Al, and ^^F{f3+f^Si, followed by 27,29,30p(^+)27,29,30gj^ One of 
the most relevant outcomes is a substantial increase in both 26'2'i'gi which, except for 28gi^ 
become the most abundant silicon isotopes in the envelope now. In fact ^^Si attains its max- 
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imum value at this stage (soon after, as the temperature rises towards Tpeafe, it begins to be 
destroyed by proton captures). Concerning ^®Si, it has significantly increased in the burning 
shell, up to 2.5 x 10~^, by mass. On the contrary, ^^Si has been slightly destroyed at the burn- 
ing shell by (p,7) reactions, whereas ^°Si remains mostly unafi^ccted. Major activity involves 
also several phosphorus isotopes: whereas ^^P remains close to its initial abundance, major 
synthesis of ^^P (reaching 4.4 x 10~^ by mass, at the burning shell) and 27,29,30p 10^^, 
by mass) is also reported. Only one sulfur isotope, ^^S, shows a noticeable variation: it is, 
in fact, increased by (p,7) reactions on ^^P nuclei at the burning shell, and shortly after, 
carried away to the outer envelope by means of convective transport. No significant nuclear 
path involving CI, Ar, K or Ca is found yet. 

When T^s reaches 3 x 10^ K (Figs. 2 to 6, fourth panel), the major nuclear activity 
is dominated by 28Si(p,7)29p(p,7)30S, ''^Si{-p,^y^P{P+)''^Si, 25Al(p,7)26Si(p,7)27p(/3+)27si, 
^^"^Al(p,7)^^Si plus ^°P(p,7)^^S. One of the most relevant outcomes related to the synthesis 
of S to Ca isotopes is the key role played by ^°P in the evolution: in fact, ^°P(p,7)"^^S 
controls the path towards heavier species, partially because the path through ^^P (i.e., 
^^P(p,7)^°S(p,7)^^Cl) is inhibited by the efficient inverse photodisintegration reaction ^^Cl(7, 
p)^°S, which is rapidly followed by a /5^-decay, leading to ^°P. Hence ^°P is a mandatory 
passing point to via 30p(p,7)3iS(p,7)32Cl(/5+)32S or through ^^P{p,^f^S{(3+f^P{p,'yf^S, 
while =*°P(/?+)3°Si(p,7)''^^P(p,7)32S is strongly inhibited by the slow ^°P /5+-decay 

Other important reactions at this stage include 2'^'269'27Al(p,7)25.27,28gj^ '^^SiiP+y^'^Al, 
28p(p,7)29S(/3+)29p(/?+)29Si(p, -ff^F and 30'3is(/5+)30,3ip^ ^ minor role is also played by 
25Si(/3+)25Al, 30p(/3+)30Si, 30S(p,7)3iCl and its reverse reaction 3iCl(7,p)30S, and both (p,7) 
and {p,a) reactions on ^^P. The outermost envelope, much cooler, is dominated by /S^-decays 
from ^^'^^Si, 29,30p marginally from ^^P and ^^K. The number of nuclei involved in the 
different dominant paths is, at this stage, very large: the most abundant silicon isotope, ^^Si, 
has already reached 3.3 x 10~^ in the inner envelope. ^^Si attains now its maximum value 
in the envelope. Regarding ^''Si, it shows a peculiar profile, with a 'hump' at intermediate 
shells. This reflects the fact that, whereas it is destroyed at the burning shell (where (p,7) 
reactions dominate synthesis from ^^P(/3+)) and at the outer part of the envelope (through 
/^"'"-decays) , it is enhanced by ^^^Al(p,7) at such intermediate shells. A similar trend is found 
for ^^Si: destroyed by (p,7) reactions in the inner envelope, and synthesized by ^^P(/5"'") at 
intermediate shells. ^'^Si increases only by a small amount. 

Concerning phosphorus, all isotopes are enhanced at this stage: 27,28,30p ^gg^ch ~ 10^^, 
by mass, whereas ^^P, the most abundant phosphorus isotope so far, achieves ~ 10^^. '^^P 
reaches ~ 10~^ at the burning shell. Except for ^'''Cl (slightly destroyed at the burning shell), 
all other chlorine isotopes are enhanced. The most abundant one, ^^Cl, reaches a mass- 
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fraction of 5.6 x 10"^. The profile of ^^Cl shows also a hump at intermediate shells, which 
distinguishes the region where destruction by ^^Cl(p,7) dominates synthesis by ^°S(p,7)^^Cl 
(innermost envelope) from that where synthesis by ^'^S(p,7)^^Cl dominates destruction by 
both proton captures and /^"'"-decays (the evolution of ^^Cl in the outermost envelope is 
dominated by ^^Cl{p+)). 

Evolution of sulfur isotopes mainly affects ^^^^^S, which show a dramatic increase pow- 
ered by (p,7) reactions on 28-30p^ jj^ fact, becomes the most abundant sulfur isotope at 
this stage. Moderate enhancement of ^^S and destruction of 33,34,36g jg g^jg^ j^cportcd. Some 
marginal activity involves also Ar and K, mainly affecting 34,35,37^^,^ enhanced at the burning 
shell, and ^^Ar, which moderately decreases. 

The course of the thermonuclear runaway drives the burning shell up to a peak tem- 
perature of Tmax = 3.26 X 10^ K (Figs. 2 to 6, panel fifth). The main nuclear activ- 
ity is now driven by ^^Al(p,7)26Si(p,7)27p and the chain ^'^Si{p,'yy^F{j3^y^Si{p,'jy^F{p, 
7)30S(/3+)30p(p,7)3iS(/5+)3ip(p, 7)32s. Significant activity also involves 24Al(p,7)25Si(/3+)25Al, 

26Si(/3+)26™Al(p,7)27Si, 269Al(p,7)27Si, ^^PiP+f'Si, 28p(p,7)29S(/3+)29p(/3+)29Si(p, 7)3°?, 

^^P(p,Q;)^^Si and ^°S(p,7)^-^Cl and its reverse photodisintegration reaction. At this stage, 
protons are energetic enough to overcome the Coulomb potential barrier of many nuclei, in 
particular ^^Si which, for the first time, begins to decrease at the burning shell. Reduction of 
2^'2^Si and a noticeable enhancement of ^^''^^Si is also found. Concerning phosphorus, 29,30p 
but mainly ^^P increase at this stage. /9+-dccays arc responsible for a significant reduction 
of both 27,28p^ ^YY^^ most relevant feature regarding the evolution of sulfur isotopes is an 
important increase of ^°~^^S. In fact, at this stage, ^^S becomes the most abundant sulfur 
isotope in the envelope. ^^Cl is favored also by the rise in temperature, and becomes now 
the most abundant chlorine isotope in the vicinity of the burning shell. 32,34g,34mQ| ^^.^ 
enhanced during this stage. Nuclear activity affects now the edge of the network, with some 
enhancement of heavier species, such as 34,35,37^^.^ 37,38j^ ^^^^^ some extent, ^^Ca. How- 
ever, the moderately low nuclear reaction fiuxes in this region justifies our choice of ^°Ca as 
the last isotope of the network. As earlier in the evolution of the nova outburst, the outer 
envelope is again dominated by /3+-decays such as 2^'2'^Si(/9"'")^^"*'^'''Al and 2^'^°P(/?+)2^'^°Si. 

Due to the sudden release of energy from the short-lived /?+-unstable nuclei ^^N, ^^'^^Q 
and ^^F, the envelope begins to expand. A hundred seconds after Tpeak, the envelope reaches 
a size of 10^ cm (Figs. 2 to 6, sixth panel). At this stage, the burning shell is dominated by 
265Al(p,7)27Si(p, 7)28p(/5+)28Si(p,7)29p, 29Si(p,7)30p(p,7)3iS, and 3ip(p,7)32s, plus a signif- 
icant contribution from 25Al(p,7)26Si(/3+)26"Al(p,7)27Si, 27Si(/5+)27Al(p,7)28Si, ^^P{(3+f^Si, 
29p(p,7)30S(/5+)30p(/3+)30Si(p, 7)3ip(p,a)28Si, and =^2,34s(p^^)33,35Q_ ^he ex- 

pansion of the envelope, which is accompanied by a drop in temperature, provides a quite 
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different chemical pattern above the burning shell. In particular, the evolution of the Si-Ca 
nuclei is fully dominated by (3~^-decays (the most important ones being ^®~^°P(/9"'")^^~^°Si, 
^'''Si(/3"'")^^A1, and ^^S(/9^)^^P, at the inner shells, whereas the outer envelope is dominated 
by 29.30p(^+)29,30gj^ 34mQi(^+)34g^ 38 j^(^+)38 . At this Stage of the evolution, the 
envelope shows dramatic changes in composition. In particular, isotopes such as 28,29,30gj^ 
30,3ip^ 32-34g^ 35q|^ 36,37^j, qj, 38j^ show a significant enhancement whereas ^^Si, ^^~^^S, 27-29p^ 
3i-34q|^ or ^^Ar are efficiently destroyed. 

During the final stages of the outburst (Figs. 2 to 6, seventh panel), convection recedes 
from the surface and the burning shell becomes detached from the major part of the enve- 
lope, which is ultimately ejected. When the envelope reaches a size of 10^'' cm, the main 
nuclear path in the Si-Ca region is basically carried by 28-30p(^+)28-30gj^ 26,27gj(^+)26m,27Ai^ 

=^°'=^^S(/3+)^°'^^P, and ^^C\((3+)^% with ^^P((3+)^^Si being the dominant one within most of 
the envelope (again, contribution from ^^"^C1(/3''")^^S and ^^K{P'^Y^Ar are important in the 
outer envelope). 

Soon after, when the envelope reaches a size of 10^^ cm (Figs. 2 to 6, eighth panel), most 
of the /^''"-unstable nuclei have already decayed. Only the medium- lived ones (r ~ minutes), 
contribute yet to the energy production: 30p(/3+)30si, 34mc^^+)34g^ ^^Kl(3+)^^Ar. The 
mean composition of the cjecta in the Si-Ca mass region consists, for this particular model, 
mainly of ^^^^''Si, ^^S, and "^^P. The most overproduced species with respect to solar arc '^^P 
and '^°Si, with overproduction factors (i.e., ratios between mean mass fractions in the ejecta 
over solar values) 1100 and 600, respectively. 32,33g ^^^^ overproduced by a factor ~ 100. 

In summary, nuclear activity in the Si-Ca mass region is powered by leakage from 
the NeNa-MgAl region, where the activity is confined during the early stages of the our- 
burst. The main nuclear reaction that drives nuclear activity towards heavier species be- 
yond sulfur is ^°P(p,7), through two different paths: either ^°P(p,7)^^S(p,7)^^Cl(/9+)^^S 
or ^°P(p,7)^^S(/3+)^^P(p,7)"^^S. Then, the nuclear flow is dramatically reduced by the slow 
^^S(p,7) limiting the production of heavier isotopes in the S-Ca mass range. 



3. Nuclear reaction rates in the Si-Ca mass region 

3.1. Proton captures on silicon isotopes 

The rates of the ^^'^^Si(p,7) reactions have been discussed in Jose et al. (1999). These 
reactions do not play a dominant role in the path towards heavier species as photo disinte- 
gration of ^^P is effective (low Q- value) and because the 27,28p life^ij^es are small (< Is). 
The ^^Si(p,7)^^P reaction is the only significant leak from the Ne-Al region towards heavier 
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nuclei (Figure 1). Hence it is natural to begin our analysis at this point, referring for the 
discussion of rates in the Ne-Al chains to our previous paper (Jose et al. 1999). Except for 
the very recent work of Iliadis et al. (2001), published after this work was completed, the two 
last compilations of thermonuclear rates end at silicon, either with ^'^Si(p,7)^^P (Caughlan 
& Fowler 1988, hereafter CF88) or 28Si(p,7)29p (Angulo et al. 1999, hereafter NACRE). 
The ^^Si(p,7)^^P rate from NACRE differs significantly from the CF88 one but only below 
10^ K, while the ^^'^°Si(p,7) reaction rates should not deviate much from the CF88 ones at 
nova temperatures (Angulo, 2000). Accordingly, in this paper we have used the ^^~^°Si(p,7) 
rates from CF88. 



3.2. Proton captures on phosphorus isotopes 

The ^^P(p,7)'^°S rate was calculated by Wiescher & Gorres (1988) from the proper- 
ties of the T=l isospin triplet. The positions of the first levels above threshold in ^'^S 
were calculated by the Thomas-Ehrman shift method while the spectroscopic factors were 
taken from the analog levels in ^'^Si. This rate is for the moment rather poorly known 
at nova temperatures as a consequence, in particular, of the uncertainty on the calcu- 
lated position of the first resonance. However, this reaction rate is not crucial as there 
exists an alternative link to heavier isotopes initiated by the ^^P /3+-dccay (Figure 1). For 
the 3°P(p,7)3iS rate (6.0371D+01, 8.6903D-01,-5.3250D+01, -4.5676D+00, 1.7458D+00,- 
1.5449D-01,-1.1437D-F01, in REACLIB format), we use the result of a statistical model cal- 
culation (Thielemann et al. 1987, hereafter SMOKER), because of the lack of spectroscopic 
information on ^^S. Indeed, even though 17 levels are known up to 1 MeV above threshold, 
very few spins are reported. Hence, it is not possible to calculate a realistic reaction rate 
based neither on ^^P analog levels nor shell model calculations (Endt 1998). However, the 
statistical model assumes a high level density which is a priori favored by the high Q-value 
(6.133 MeV) but remains questionable at the moderate temperatures found in novae. Hence, 
with such a level density, it is more likely that the rate is governed by the properties of a 
few individual levels. 

The inspection of Figure 1, shows that this reaction plays a key role as ^''P is a mandatory 
passing point to reach heavier nuclei and with its relatively long lifetime (ti/2=2.5 mn) it stops 
further nucleosynthesis unless proton captures arc fast enough. Accordingly we investigate 
the effect of a possible change of this rate by a test factor of 100 (sec Section 4). The 
^'^P(p,7)^^S and ^^P(p,Q;)^^Si reaction rates are deduced from experimental results (Iliadis et 
al. 1991; 1993; Ross et al. 1995). While the first one is well known, the second still suffers 
from uncertainties (^100) between Tg — 0.1 and 0.2 due to unknown resonance strengths. 
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However, even the upper limit of the ^^P(p,Q;)^^Si rate remains smaller than the ^^P(p,7)^^S 
one at nova temperature. Hence this uncertainty should scarcely affect ^^P destruction and 
nucleosynthesis beyond phosphorus. 

3.3. Proton captures on sulfur isotopes and beyond 

For ^°S(p,7)^^Cl, we use the rate calculated by Herndl et al. (1995) from shell model 
calculations. Another rate was provided by Van Wormer et al. (1994) based on the properties 
of the analog levels in ^^Cl. The two rates differ by a factor of ^30 (Fig. 1 in Herndl et 
al. (1995)) due to the unknown location of the first excited state in ^^Cl. This should 
however have little effect since photodisintegration of ^^Cl is fast (Q=0.286 MeV). The 
^"^S(p,7)^^Cl reaction rate was first studied by Vouzoukas et al. (1994). They measured 
the energy of the relevant ^^Cl levels via a ^^S('^He,t)'^^Cl transfer reaction. They used 
experimental spectroscopic factors from the ^^P analog levels and radiative widths from shell 
model calculations and their calculated ratios (F^/r^) are compatible with the Lefevbre et al. 
(1997) measurements. A detailed investigation of the uncertainty affecting this rate has been 
performed recently by Ihadis et al. (1999) leading to a nominal rate (adopted here) slightly 
lower than the previous one (Vouzoukas et al. 1994), and a factor of 3 uncertainty in the 
domain of nova nucleosynthesis. One can thus consider this rate as sufficiently known. The 
^^S(p,7)^^Cl reaction rate was calculated by Iliadis et al. (1992) from measured resonance 
strengths or extracted spectroscopic factors (see also Iliadis et al. 1999, footnote in Table 
7). The uncertainty on this rate was found to be very small (Thompson & Iliadis 1999). It 
is important to notice that the rate deduced from experimental data is much smaller (see 
Ihadis et al. 1992) than the previously adopted reaction rate. 

For the reactions ^^S(p,7)^^^'^^"^Cl, we have taken into account the formation of both 
the ^^Cl ground (ti/2 = 1.53 s) and isomeric {ti/2 — 32.0 min) states which, as in ^^Al, have 
to be treated as separate isotopes (Coc et al. 2000a) at nova temperatures. We considered 
the resonances corresponding to the first five levels above threshold (Endt 1990). Their 
branching ratios to the ground or isomeric states were extracted from Endt (1990). For the 
four unknown resonance strengths, we used upper limits with the usual 0.1 reduction factor. 
They were obtained from estimates of the proton widths (Fp ^ F^) based on calculation of 
single particle in a Wood-Saxon potential whose depth is adjusted so as to reproduce the 
excitation energy of the resonance. The corresponding estimated ^^S(p,7)^"^™Cl rate is given 
by: 

NA{crv) = 1.5923x10^ T~^^^ x (0 to 1) x 

(2.4x10-29 exp (-0.3389 /Tg) x 1. 
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+1.5x10-1^ exp (-2.003 /Tg) x 0.89 
+2.4x10-5 exp (-2.836 /Tg) x 0.64 
+1.X10-3 exp (-4.623 /Tg) x 0.76) 
+2.13x10^ Tg-o-^^'^^exp (-4.814 /Tg) x 0.5 

where the multiplying factors (1., 0.89, 0.64, 0.76 and 0.5) correspond to the branching 
ratios to the isomeric state (Endt 1990). The ^'^S(p,7)^'*^Cl is hence given by replacing these 
factors by (0., 0.11, 0.36, 0.24 and 0.5). The last term is a multiresonant contribution from 
the higher energy resonances (Endt 1990) with an assumed average branching ratio of 0.5. 

Formation of ^^""Cl proceeds only through the ^^S(p,7)^^'"Cl reaction since, after internal 
transitions, "^^Ar /3+-decays lead only to ^^^Cl (a similar configuration is found for ^^Al). An 
important aspect of ^^"^Cl nucleosynthesis is its destruction by thermally induced radiative 
transitions to its short lived ground state. The corresponding effective lifetime, that we use 
here, has been calculated (Shell Model) by Coc et al. (2000a): it exhibits a rapid evolution 
(3 orders of magnitude) from its laboratory half-hfe value of 32.0 min down to ^1 s. This 
transition occurs around Tg RiO.15-0.2, a typical domain of temperatures in novae, and this 
strongly modifies the prospect for its detectability through its gamma-ray signature. Other 
nuclei, such as ^"^Al, have also isomers in the mass region of nova nucleosynthesis but they 
are at the fringe of nuclear flow and their lifetimes are much less affected than ^^Cl at nova 
temperatures. 

With some notable exceptions like ^^Cl{p,^) and ^^Ar(p,7) (see also the very recent 
compilation of Iliadis et al. 2001), experimental and theoretical (i.e.. Shell Model, see Herndl 
et al. 1995) data are in general not sufficient to allow the explicit calculation of the rates 
at nova temperatures for reactions beyond sulfur. Accordingly, we use the SMOKER rates 
instead (Thielemann et al. 1987) which are virtually identical to the more recent MOST 
rates from Goriely (1996), for these nuclei close to stability. 

4. Results and discussion 

Table 2 shows the mean composition of the cjccta, for the Si-Ca mass region, correspond- 
ing to the 1.35 Mq ONe nova described in Section 2 (hereafter. Model 135nom). The most 
abundant nuclei, for this specific mass range, are ^^S, ^^Si, ^°Si, ^^P and ^^Si. For the purpose 
of comparison, a very similar model (i.e.. Model ONel35B. See Jose, Coc & Hernanz 1999), 
computed with a previous nuclear reaction network, is also shown in Table 2. Both nuclear 
reaction networks are very similar up to silicon. Recent updates included in the new com- 
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putations involve ^^F(p,7), ^^F(p,ci;) (Coc et al. 2000b, nominal rates), ■^'''0(p,7), ^'''0(p,Q;) 
(NACRE, nominal rates), and also ^^P(p,7), 3^P(p,Q;), ^^S(p,7), ^^S(p,7), ^^S(p,7) (leading 
to both 345,34m(2;j^_ ^^^i ^]^g corresponding 34g,34mQj^ temperature-dependent /^"'"-decays (see 
Section 3 for specific references). As expected, no significant differences are found for Si or 
P yields when we switch to the updated network since, as discussed in Coc et al. (2000b), 
proton captures onto ^^O or ^^F don't play a dominant role in the energetics of the explo- 
sion, and hence, a similar peak temperature and characteristic timescales for the evolution 
are found. Differences, however, are found for ^^S and beyond. In fact, the new ^^S yield 
increases by a factor of 2, because of the larger new ^^P(p,7) rate (below 10^ K and also in 
the vicinity of Tpeak) and the lower ^^S(p,7)^^Cl destruction rate. On the contrary, ^^'^^S, 
35,37q| g^j^^ 36,38 dramatically reduced by, at least, a factor of 10. This remarkable 
change results mainly from the update of the ^^S(p,7) rate, much lower than the previous 
one in the domain of nova temperatures, as pointed out by lliadis et al. (1992). Hence, 
whereas the arguments to invoke a likely nova origin for a number of presolar mctcoritic 
grains (sec Amari et al. 2001) remain unaffected by the results reported in this paper (i.e., 
sihcon excesses), the differences found for other nuclear species, in particular the ^^S/^^S 
ratios, are relevant for the unambiguous identification of future nova grain candidates: for 
instance, the previous ^^S/^^S isotopic ratio of 3.1, found for Model ONel35B, has now raised 
up to 97 for Model 135nom, basically due to the update of the ^^P(p,7) and ^^S(p,7) rates. 

It is also worth noting that our numerical calculations reveal the presence of the gamma- 
ray emitter ^'^'^Cl in the ejected envelope. However, its content has been dramatically reduced 
by the update of the =^2s(p,7)33Cl and 33s(p,7)34m,34gci ^^^^^ together with the =^^"*C1 
decay: whereas a mean mass fraction of 3.4x 10~^, 30 minutes after Tpeak, was found for Model 
ONcl35B, only 7 x 10^^ by mass, remain, for the same time, in the present Model 135nom. 
Hence, the possibility of a potential detection of the gamma ray emission corresponding to 
34mQj f3~^-decaj, invoked as another signature of a classical novae (see Leising & Clayton 
1987), seems to be ruled out with the present calculations. 

As pointed out, one of the most relevant outcomes of this study, from the nuclear physics 
viewpoint, is the crucial role played by '^°P(p,7) in the synthesis of heavier nuclear species. 
Its purely theoretical rate is questionable since the model used for its calculation reaches 
its limits for such a light nucleus at low temperature. Hence, we have estimated its impact 
on the yields by computing two new evolutionary sequences of nova models, for which we 
have arbitrarily modified the ^°P(p,7) rate: Model p301ow, where the ^°P(p,7) has been 
reduced by a test factor of 100 throughout the calculation, and Model p30high, where the 
nominal rate has been arbitrarily multiplied by 100. Results are summarized also in Table 
2. The increase by a factor of 100 affects only the ^°Si yield: since the path towards ^^S is 
favored (which in turn, accounts for the moderate increase of ^^P in the ejecta), much less 
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^°Si (coming from ^°P(/3+)) is left in the envelope. However, the impact on other nuclear 
species is very limited. On the contrary, when the ^'^P(p,7) rate is reduced by a factor of 
100, dramatic changes in the expected yields occur. Now, ^°P(/3+) competes favorably with 
proton captures, and an important enhancement of '^"^Si is obtained. In fact, the reduction 
of the proton capture rate halts significantly the synthesis of elements above Si (by a factor 
of ~ 10 with respect to the values found with the nominal rate). We have compared our 
results with the latest available results from hydrodynamic simulations of ONe novae by 
Starrfield et al. (2001): they state that the switch to an updated nuclear network (Iliadis 
et al. 2001) is accompanied by an important reduction of the abundances of nuclei above 
aluminum, in good agreement with the results found here. However, our new computations 
show that reduction of the yields affect only nuclei above '^^S. In fact, a significant reduction 
of ^^S results only when we adopt the low ^°P(p,7) test rate. Hence, this discrepancy might 
be likely attributed to a different prescription for the adopted ^°P(p,7) rate. In particular, 
the rate used by Starrfield et al. (2001) seems to be somewhat lower than the nominal rate 
adopted in our paper, which stresses again the crucial dependence of the S-Ca yields on this 
particular rate. 

Despite a detailed comparison with observations is out of the scope of this paper, because 
a larger number of nova models is required to properly explore the wide parameter space, 
some conclusions can be outlined. In fact, spectroscopic abundance determinations show 
evidence for some overproduction of nuclei in the Si-Ca mass region in a number of classical 
novae, including silicon (Nova Aql 1982, Snijders et al. 1987, Andrea et al. 1994; QU Vul 
1984, Andrea et al. 1994), sulfur (Nova Aql 1982, Snijders et al. 1987, Andrea et al. 1994), 
chlorine (Nova GQ Mus 1983, Morissct & Pequignot 1996), argon and calcium (Nova GQ 
Mus 1983, Morisset & Pequignot 1996; Nova V2214 Oph 1988, Nova V977 Sco 1989 and 
Nova V443 Set 1989, Andrea et al. 1994). Indeed, models of explosions in 1.35 M© ONe 
white dwarfs have an end-point for nucleosynthesis below calcium, hence in agreement with 
observations. This suggests that peak temperatures attained at the burning shell during nova 
outbursts cannot be much larger than 3 x 10^ K, otherwise an rp- process would result at 
such high temperatures, driving the nuclear path towards heavier species, not observed so far. 
Moreover, as stated in the discussion of the main nuclear paths, determination of chemical 
abundances in the Si-Ca mass region during nova outbursts is affected by uncertainties 
of nuclear physics origin. In particular, the uncertainty associated to ^''P(p,7) makes it 
difficult to determine realistic abundances above silicon, and to clarify as well whether the 
path towards heavier species is halted or, on the contrary, if significant production of Si-Ca 
nuclei is still expected. A better determination of this crucial rate would be of significant 
importance for the studies of nova nucleosynthesis that takes place in massive ONe white 
dwarfs. 
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Fig. 1. — Main nuclear paths in the Si-Ca mass region. 
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Fig. 2. — Snapshots of the evolution of sihcon isotopes (in mass fractions) along the accreted 
envelope, for a 1.35 Mq ONe nova accreting at a rate of M = 2 x 10^^° MQ.yr"^ . The mass 
coordinate represents the mass below the surface relative to the total mass. From top to 
bottom, panels correspond to the time for which the temperature at the burning shell reaches 
8 X 10'^, 10^, 2 X 10^, 3 X 10^, and Tpeak = 3.26 x 10^ K, plus three panels corresponding to 
the last phases of the explosion, when the white dwarf envelope has already expanded to a 
size of Rwd ~ 10^, 10^°, and 10^^ cm, respectively. The base of the ejected shells correspond 
to a mass coordinate of —3.05 x 10~^. 
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Fig. 3. — Same as Fig. 2, for phosphorus isotopes. 
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Fig. 4. — Same as Fig. 2, for sulfur isotopes. 
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Table 1: Initial composition (Si to Ca) in the accreted envelope, assuming a 50% degree of 
mixing with the ONe white dwarf core (see Ritossa et al. 1996). 



Nuclei 


Mass fraction 


28si 


3.3E-4 


29Si 


1.7E-5 


30Si 


1.2E-5 


31p 


4.1E-6 




2.0E-4 


33S 


1.6E-6 


34g 


9.3E-6 


36S 


4.7E-8 


35Q 


1.3E-6 


3^C1 


4.3E-7 


36 Ar 


3.9E-5 


38 Ar 


7.7E-6 


39K 


1.7E-6 
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Table 2: Mean composition of the ejecta (Si to Ca) from 1.35 Mq ONe white dwarf models. 
Solar values are given for comparison. 



Nuclei 


Solar 


Model 


Model 


Model 


Model 






ONcl35B 


p301ow 


135nom 


p30high 


28gi 


6.5E-4 


4.2E-2 


3.9E-2 


4.4E-2 


4.5E-2 


29Si 


3.5E-5 


2.1E-3 


1.9E-3 


1.9E-3 


1.9E-3 


3°Si 


2.4E-5 


1.4E-2 


6.8E-2 


1.4E-2 


4.6E-4 


31p 


8.2E-6 


8.4E-3 


l.lE-3 


9.3E-3 


1.2E-2 


32g 


4.0E-4 


2.3E-2 


4.6E-3 


4.5E-2 


5.5E-2 


33S 


3.3E-6 


7.6E-3 


4.6E-5 


4.8E-4 


5.9E-4 




1.9E-5 


9.3E-3 


2.2E-5 


2.2E-4 


2.7E-4 


36S 


6.4E-8 


6.6E-9 


5.6E-9 


6.0E-9 


6.1E-9 


35Q 


3.5E-6 


4.5E-3 


2.9E-5 


2.9E-4 


3.6E-4 


37C1 


1.2E-6 


1.4E-4 


9.0E-6 


1.7E-5 


2.0E-5 


36 


7.7E-5 


2.2E-3 


5.1E-5 


1.9E-4 


2.3E-4 


38 A^ 


1.5E-5 


2.7E-5 


8.1E-6 


9.1E-6 


9.5E-6 


39K 


3.5E-6 


5.2E-6 


4.8E-6 


4.7E-6 


4.7E-6 



